Table 1 

Values of the hadronic masses, measured with 
kvaience = k^ea = 0.1570, usiug APBC ou Spa- 
tial directions on the sea quarks, PBC on the 
valence quarks and a 7^ smeared source. We 
report also the value for the square ratio of the 
pseudoscalar to the vector mesons and the value 
for the plaquette. 
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ferent methods: A) we study the behavior of the 
jacknife error a^ (calculated by dividing the mea- 
sures in N bins of size T) as a function of the 
size T, analyzing when the error reaches a plateau 
and, B), we calculate the autocorrelation function 
C(t) at time r, (see eq.(2.19) of ref.[7]). 

With both methods we obtain for the plaquet- 
te and the hadrons correlation functions, at Eu- 
clidean time 16 (the middle of the temporal lat- 
tice), an autocorrelation time of about 50 HMC 
trajectories (see fig. 2). Starting from these re- 
sult we estimate our statistical errors with the 
jacknife method used on 10 bin of 5 measures (50 
trajectories) each. Results are reported in Table 
1. We can see from this table that the masses are 
estimated with a statistical error of about 1%. 

In order to study the systematic errors coming 
from the finite lattice, we calculated the hadronic 
masses with both PBC and APBC on valence 
quark. The difference between the values of 
masses obtained gives us an estimation of the con- 
tribution from < P >^ of eq. 1. From fig. 3, 
we can see an error of 3% on the estimation of the 
plateau of the tt and of 6% for the p. The sys- 
tematic finite size errors are greater than the sta- 
tistical errors quoted in Table 1. This difference 
is lower than that found for staggered fermions 
in fig 2 of ref.[2] as expected since our lattice 
is almost a factor 2 larger. However the effect 
is in the same direction: if the boundary condi- 
tions for sea quarks are AP, the masses obtained 
with PBC on valence quarks are lower than those 
obtained with APBC on valence quarks, as ex- 



Figure 3. Pseudoscalar effective mass with both 
PBC and APBC on valence quarks, from prop- 
agators without smeared source. 



pected from eq. (1). This is true for both Wilson 
and staggered fermions. 

We want to stress that the smearing procedure 
is crucial also in full QCD, as we have observed 
in the study of the baryons correlation functions. 
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ence of the fermionic part of the QCD action is 
analogous to the existence of a magnetic field h 
which breaks the Zs symmetry in a model for a 
spin, n, which can take the three possible Zs val- 
ues, n is coupled to the external magnetic field 
via the Hamiltonian 



Hh = hU + h^U^ 



(2) 



Hh is not Zs invariant. The values taken by H 
are in correspondence with the expected phases of 
Polyakov loops, while the values of h are in cor- 
respondence with the chosen sea quarks bound- 
ary conditions. Following this model we expect 
that with PBC on sea quarks the e'^'^'^ and 
the e~'^'^'^ phases of the Polyakov loops will be 
present with equal probability, while Polyakov 
loops are likely to be polarized in the Hq direction 
if APBC are used. Similarly we expect to be able 
to align the Polyakov loops along the e'^'^'^ (or 
g-827r/3-j -j^ ^j^g ^■g space, if we choose —e-'^f/s 
(or _e'2T/3 respectively) boundary conditions on 
the sea quarks. 

We checked this suggestions with a simula- 
tion, performed on APEIOO, with 2 fiavor Wil- 
son fermions at /3 = 5.3 on a 8^ x 32 lattice with 
kgga = 0.1670. We carried out two different runs, 
one with fully PBC on the sea quarks and the 
other one with APBC in the spatial directions 
and PBC in the temporal one. Gauge configura- 
tions have been produced with APEIOO using the 
HMC algorithm^. We created a set of 440 (ther- 
malization) -|- 1350 thermalized trajectories. We 
have performed the measurement of the spatial 
loops, Px,Py and P^,, on one every 5 consecutive 
trajectories, using the smearing procedure^®] to 
reduce fiuctuations. With APBC on spatial di- 
rections, the phases are close to zero, while with 
PBC, the phases are concentrated in two regions, 
near e'^'^'^ and g"'^'^'^, as we can see from fig. 1. 

Fixing the Polyakov loop phase reduces the 
statistical fiuctuations on the hadron propagator 
and, hence, on the computed hadron masses. Fur- 
thermore the use of PBC on sea quarks causes 
a flip-flop effect between the two phases of the 
Polyakov loops, g"'^'^'^, g'^'^'^, which can pro- 
duce a longer thermalization time. 



2. SIMULATIONS ON 16^x32 AT /3 = 5.55. 

We started with a run on a 16^ x 32 lattice 
volume at /3 = 5.55 with kg^a = 0.1570. We 
use APBC on spatial directions on sea quarks. 
We created a set of 900 (thermalization) -|- 500 
trajectories. The last 500 have dt = 0.083 and 
Nmd = 60 (dtxNMD ~ 0.5). We have performed 
a measure, every 10 trajectories, of the hadronic 
correlation functions, of the 0"*""*" and 2+"'" glueball 
correlation functions and of the Polyakov loop. 
For each trajectory we estimate the average value 
of the Wilson loop, Wixi, (the plaquette). In the 




Figure 2. Autocorrelation for the pseudoscalar, 
TT, meson correlation functions at the Euclidean 
time t = 16. On the top we present the results 
from method (A) and on the bottom the ones 
from method (B). 



calculation of the valence quark propagator we 
use a smeared source of size 7^, while for the glue- 
ball correlation functions we use 21 iterations of 
the smearing. The physical lattice size (estimated 
from ref.[6]) is _L ~ 1.5 fm. 

The first goal that we want to reach is an un- 
derstanding of the correlations between measures 
done on consecutive HMC trajectories. This stu- 
dy is important in order to estimate a correct 
value for the statistical error. We use two dif- 
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1. THE SPATIAL POLYAKOV LOOPS 
PHASE. 
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Polyakov loops are responsible for the differ- 
^^ ence between quenched and unquenched finite 

size effects on the QCD mass spectruml^lPl. The 

reason can be understood by looking at the va- 
T-H lence quark hopping parameter expansion of the < 

^T meson propagator [2]: 

0\ 

"^ G = ^fcl^)MC)> + ^fcl^V,„,(P(C)> . (1) 

I 20 1 

rr^ The sums extend over all closed paths (C) of 'X'*'^^ 

(— I length 1(C). W(C) are standard Wilson loops, lo r 

completely contained into the lattice, while P(C') 5 r 

are Polyakov type loops and (•) denotes gauge 

field average; the value of the index ayai depends „ 

on the spatial boundary conditions on the va- 
lence quarks: ayai = +1 for periodic boundary 
conditions, PBC, and ayai = (~1)" for antiperi- 
odic boundary conditions, APBC, respectively, 

with n the number of windings around the lat- Figure 1. Histogram of the phase, 9, of the x 

^ji,g component of the Polyakov loop, P^^ , for APBC 

In full QCD the fermionic part of the action and P^ , for PBC on sea quarks, 

(both Wilson and staggered fermions) breaks ex- 
plicitly the simmetry that, in quenched confined 
QCD guarantees < P >= 0. A simple modeU^^ 
may be useful to clarify the situation. The pres- 
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